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Super-hydrophobic array of grooves containing trapped gas (stripes) , have the potential 
to greatly reduce drag and enhance mixing phenomena in microfluidic devices. Recent 
work has focused on idealized cases of stick-perfect slip stripes. Here, we analyze the 
experimentally more relevant situation of a pressure-driven flow past striped slip-stick 
surfaces with arbitrary local slip at the gas sectors. We derive approximate formulas for 
maximal (longitudinal) and minimal (transverse) directional effective slip lengths, that 
are in a good agreement with the exact numerical solution for any surface slip fraction. 
By representing eigenvalues of the slip length-tensor, they allow us to obtain the effective 
slip for any orientation of stripes with respect to the mean flow. Our results imply that 
flow past stripes is controlled by the ratio of the local slip length to texture size. In case of 
a large (compared to the texture period) slip at the gas areas, surface anisotropy leads to 
a tensorial effective slip, by attaining the values predicted earlier for a perfect local slip. 
Both effective slip lengths and anisotropy of the flow decrease when local slip becomes of 
the order of texture period. In the case of small slip, we predict simple surface-averaged, 
isotropic flows (independent of orientation). 



1. Introduction 

The develop ment of microfluidi c s has motivated interes t in manipulating flows in very 
small channels ( Stone et al. 2004 ; Squires fc Quakd 2005 ). Most of microfluidic devices 
operates with a pressure flow, which is faced with two main difficulties at this scale and 
under typical operating conditions. First, it is difficult to drive such a flow due to huge 
hydrodynamic resistance. Second, it is very difficult to mix, which normally requires a 
generation of a tranverse flow. 

An efficient strategy for moving efficiently fluid in a tiny channel is to exploit hy- 
drodynamic slip, which can be generated at hydrophobic surfaces and is quantified 
by the slip leng th b (the distance within the solid a t which the flow profile extrap- 
olates to zero) ( Vinogradova 19991 Lauga et all 12007 ; Bocauet fc Barrat] 2007 ). Since 
for hydroph obic smooth and homogeneous surfaces b can be of the order of tens of 
nanometers (Vinogradova fc Yakub"ovll2003t Cottin-Bizonne et al. 2005 : Jolv et al. 2006j 
Vinogradova et aZ.ll2009h . but not much more, it is impossible to benefit of such a nano- 
metric slip for pressure-driven microfluidic applications. However, super-hyd rophobic 



(SH) textures can significantly am plify hydrodynamic slip due to gas entrapment ([Vinogradova et al 
19951 : Cottin-Bizonne et al. 2003f ) leading to the huge slip length at the gas area. The 
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Figure 1. (Left) Sketch of SH stripes: 6 = n/2 corresponds to transverse, whereas 9 — to 
longitudinal stripes; (right) situation in (left) is approximated by a periodic cell of size L, with 
equivalent flow boundary conditions on gas-liquid and solid-liquid interface. 



composite nature of the texture, however, requires regions of lower slip (or no slip) in 
direct contact with the liquid, so the effective slip length of the surface b c g is reduced. 
Indeed, experimental stu dies of flow past SH surfaces suggest that effective slip is of the 
order of several microns ( Ou fc Rothstein 2005 : Joseph et aLll200rj IChoi et aZll2006f l. 

SH surfaces consisting of periodic array of grooves containing trapped gas (Cassie's 
state) are especially interesting since they allow to highlight effects of anisotropy. For 
anisotropic textures 6 e ff varies with the orientation of the wall texture relative to flow 
and is generally a tensor ( Bazant fc Vinogradovall2008l) . Such surfa ces have been already 
used for reduction in p ressure-driven flows ( Ou fc Rothstein! 2005 ) and enhancement of 
mixing ( Ou et all 12007 ). The problem of flow past stripes has been examined theoreti- 
cally mostly with a typical geometry sketched in Fig. [T] corresponding to of a roughly 
flat (no meniscus curvature) liquid interface, so that the modeled SH surface appeared 
as a perfectly smooth with a pattern of boundary conditions. In the case of thin chan- 
nels {H <C L, where H is the channel thickness, and L is the period of the texture) 
the problem was solved for any two-component (e.g. low-slip and high-slip) texture, and 
striped surfaces were shown to provide rigor ous upper and lower bo unds on the effec- 
tive slip over all possible two-phase patterns (jFeuillebois et al. l2009h . The quantitative 
understanding of liquid slippage past such a surface in the thick channel (H ^> L) is 
however still challenging. Press ure-driven flow has been analyzed for an idealized case of 
a perfect slip at the gas area (Philip 1972 : Lauga fc Stond 2003t Cottin-Bizonne et all 
2004 ISbragaglia fc Prosperettill2007t ) and led to 
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where 4>2 — S/L denotes the fraction of the liquid-gas interface (correspondingly, <j)\ = 
1 — 02 is the fraction of solid-gas area), with the typical length scale of the slipping area 
(5, and b^a and fr|| ff denote ef fective transverse and longitudinal slip lengths. Following 
iBazant fc Vinogradova ( 2008 ) , these are the eigenvalues of the second-rank effective slip- 
length tensor b c s represented by symmetric, positive definite 2x2 matrix diagonalized 
by a rotation: 

b cff = S e ( ^ L ° ) S_ 9 , S e = ( t a ) (1-2) 



Therefore, Eqs. (jl.ll) allow to calculate b c g in any direction given by an angle 9 (Fig.[T]). 

Eqs. (|l.ll) provide an upper limit for the effective slip lengths and in many situations 
would be expected to overestimate the m. One reason is the possible meniscus curvature, 
which has been clarified in recent work (|Sbragaglia fc Prospere tti 2007; Hvvaluoma fc Hartina 
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2008; Davi s fc Laug a 2009). Another is viscous dissipation t aking place in the u nderlying 
gas phase. Indeed, the more realistic "gas cushion model" ( Vinogradoval 1995 ) predicts 
the finite slip length at the slipping area 



(1.3) 



where e is the thickness of the gas layer, rj is the viscosity of liquid, and r/ g is the viscosity 
of gas. Taking into account that under typical conditions p/rif, ~ 50, t he variation of the 
SH texture height, e, in the typical interval 0.1 — 10 /jm (jQuere l2005h gives b — 5 — 500 
/im, i.e. b might be as small as typical L or even less. For this reason, it is attractive 
to consider this experimentally relevant situation. However, despite its fundamental and 
practical significance, pressure-driven flow o ver partial slip stripes has received little 



attention. This has been studied numericallvf Cottin-Bizonne et a/.ll2004t iPrieziev et al 



20051: lYbert et aLll2007h . Nevertheless, no analytical resolution of the Stokes equation 
with this set of boundary conditions has been performed up to now. 

In this paper, we provide analytical solutions to pressure-driven flows over SH stripes. 
In §2 we formulate the problem and derive expressions for the effective slip for longitu- 
dinal and transverse stripes, which allows us to obtain a solution for any orientation of 
stripes with respect to a gradient of pressure. In §3 we compare our results with numerical 
calculat ions performed by C. Cottin -Bizonne and C. Barentin using the method devel- 
oped in ICottin-Bizonne et all ( 2004 ) and discuss implications for the use of SH stripes 
to control hydrodynamic flows. We conclude in §4. 



2. Model and Analysis 

We consider a pressure-driven flow past an idealized, flat, periodic, striped SH surface 
in the Cassie state (sketched in FigO}, where the liquid-solid interface has no slip (&i = 0) 
and the liquid-gas interface has partial slip (b 2 = b) . Our results apply to a single surface 
in a thick channel (H ^> max{L,6}), but not to thin channels (H <C m injL, b}) where 
the effective slip scales with the channel width ( Feuillebois et al\ \2009). The origin of 



coordinates is placed it the plane of liquid-gas interface above the middle of the slot. 
The a;-axis is defined al ong the pressure gradient, whi le the y-axis is aligned across the 



channel. According to IBazant fc Vinogradoval ([20081 ) the general problem reduces to 



computing the two eigenvalues, &jj ff an( i fyrfr> which attain the maximal and minimal 
directional slip lengths, respectively. 

The fluid flow satisfies Stokes' equations 

r)V 2 u = V P , V-u = 0, (2.1) 

where u is the velocity vector, and the applied pressure gradient is parallel to the x axis 
direction: 

Vp o = (-a,0,0) (2.2) 
The slip boundary conditions at the channel walls are defined in the usual way: 

<9u 

u(x,0,z) = b(x,z) ■ — (x,0,z), y ■ u(x, 0, z) = 0. (2.3) 
dy 

dw 

u(x, H, z) = -b H ■ — (x,H, z), y • u(x, H, z) = 0. (2.4) 
dy 

Here the local slip length b(x, z) is generally the function of both x and z coordinates. 
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For simplicity, we now consider here the case bn = 0. As the problem is linear in u, we 
seek the solution in the form: 

u = u + ui, (2.5) 
where Uo is the velocity of the flow over the homogeneous plane with the no-slip condition: 

17 



u = (^0,0,0), u 
du 
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(2.7) 



and U! is the perturbation of the flow, which is caused by the presence of the texture 
and decays far from the bottom of the channel. 

We are interested in the effective slip length b c g of the lower surface defined as 



f du\ 

\dy) ! 



(2. 



where 



means the average value in plane xOz. 



2.1. Longitudinal stripes 

In this case the problem is homogeneous in a;-direction (-Jj =0). The slip length b(x, z) = 
b{z) is periodic in z with period L. The elementary cell is determined as b(z) = b at 
\z\ < 8/2, and b(z) = at 8/2 < \z\ ^ L. In this case velocity Ui = (ui,0,0) has only 
one nonzero component, which can be determined by solving the Laplace equation with 
the boundary conditions discussed above. By choosing L/(2n) as the length scale and 
<jL 2 i '(Aj^rf) as the velocity scale we obtain in the dimensionless form 



ui(y,z) 



a n cos(nz)e 



(2.9) 



n=l 



(The sine terms vanish due to symmetry.) Condition (|2.3j) leads to the dual trigonometric 
series 
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^0 \ , v 

— + a n COS(nz) =0, C < Z %. 7T, 



(2.11) 



n=l 



where c = 7T02 and Co = Cg • 27T7// (uL) = ixH/L. To solve these series we assume that 



«o 

y 



+ a n cos(nz) = cos(z/2) / 
„=i ^ 

According ( Sneddon] 1966t ) we then get 
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Vcos z — cos f 
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< z < c. 



n = 1,2,3,. 



(2.12) 
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where P n is Legendre polynomial, and one can then show that the effective slip length 
is given by 

II ^oo ( } 

cff 2tt2Co V ; 
By integrating (|2 . 10|) in the interval [0,z], and substituting (|2.12p and (|2.14|) we obtain 

(0 < z c) 



2tt6 



h(t)dt z 

, = sec - 

V cos t — cos z 2 
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We further change the order of integration in parentheses to get 
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h(t)dt 



Vcos £ — cos t 
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The evaluation of (|2.16p gives 

r COS cZ«g 



a/cos I — COS t 



V2- 



y/cOS £ — COS t 

(2.17) 
(2.18) 



so that we get 
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h(t)dt z 
—^^ = ^^ = = sec - 

V cos t — COS 2 2 



2tt6. 
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This can be simplified by neglecting the last term in parentheses, which is small as 
compared to the main term (due to properties of arccos ^^-rj) and, thus, 
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(2.22) 



2.2. Transverse stripes 

In this case the pressure gradient depends on x, so that it is convenient to introduce 
a stream function ip(x,y) and the vorticity vector oj(x,y). The two-dimensional ve- 
locity field corresponding to the transverse configuration is represented by u(x, y) = 

yBy > — 75-c ' 0) ' an d the vorticity vector w(x, y) = V x u = (0, 0, w) has only one nonzero 



6 A. V. Belyaev and O. I. Vinogradova 




-? 6 ? -7 2 

10 10 10 10 10 10 

b/L b/L 



Figure 2. (Left) Eigenvalues fe" ff (dashed curve) and b^r s (solid curve) of the slip-length tensor 
b e fi for stick-slip stripes of period L and slipping area fraction cj>2 = 0.5 as a function of the 
local slip length b of this area. Dash-dotted curve represents the effective slip in the direction of 
driving force for tilted (8 = n/4) stripes. Symbols show numerical results. (Right) The ratio of 
theoretically predicted eigenvalues of the slip-length tensor b e a (solid curves) and corresponding 
results of numerical modeling (symbols). From left to right, <f>2 = 0.05, 0.5, and 0.95. 



component equal to 

uj = -VV (2.23) 

The solution can then be presented as the sum of the base flow with homogeneous no-slip 
condition and its perturbation due to the presence of stripes 

tp = ^ + tpi, w = fl +<•>!, (2.24) 
where and correspond to a typical Poiseuille flow 

*o = --^ + C *^, n = -y-C*. (2.25) 
r\ 6 2 rj 

The problem for perturbations if}\ and wi of the stream function and z-component of 
the vorticity vector reads 

V 2 ^i = -oji, VV = 0, (2.26) 
which can be solved by applying boundary conditions (|2.3p and (|2.4p . that take form 
( Prieziev eraZll2005l) : 

^(x,y = 0) = b(x) ■ [C * - u 1 (x,y = 0)] , (2.27) 

^(x,y = H)=0, (2.28) 
ay 

and an extra condition that reflects our definition of the stream function 

^ 1 {x,y = Q>)=Q. (2.29) 

This can be solved exactly to get 

oo 

Ui(x, y) = + ^a n cos{\ n x)e~ XnV , (2.30) 

1 n=l 

ip 1 ( X7 y) = -^ y 2 + (3 oy + Y^ (f3 n + ^-jL\ cos{XnX )e-^y, (2.31) 
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where A n = (2im)/L is the wave-number. Condition (|2.28|l leads to /3q = aoH/2, and 
(I2~29l) gives Pn=0. 

Applying boundary conditions, we obtain another dual series, similar to (12.10)) and 

(Em 



00 1 + 



H 



2tt6 



2-Kb 



a n I 1 + 2 • —j—n I cos(nx) = ——Co, < x < c, 



(2.32) 
(2.33) 

^2 — "» - 2n ffi . ( 2 - 34 ) 

and b^ s = (L/27r)(ao/Cb). Since fr/if is negligibly small, the dual series can be simplified 
to obtain 



Here 



do + /J a„ cos(na;) =0, c < x ^ 
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2.3. Arbitrary direction 

Here we consider the situation when the pressure gradient is aligned at the angle to 
the stripes. The surface velocity u s = (u s ,0,w s ) has only two nonzero components. We 
establish the coordinat e syst em so that — Vpo is parallel to the :r-axis. According to 
Bazant fc Vinogradoval (|2008l ) 



J cff • 



du 

dy 



(2.36) 



where b c g is given by Eq. (jl.2p . Average components of surface velocity then read 

-fhW 



(u.) = (bl s cos + Ks sin" 8) • C*, 



(2.37) 



(w s ) = (bl s ~ b^ s ) sine cos0 • C *. (2.38) 

The absolute value of the slip velocity on the striped SH-surface |U S | and the angle ip 
between the driving force (— Vpo) and U s are then given by 



( b ls ~ frjjf) sing cos < 
(&jj ff cos 2 + 6^ ff sin 2 i 



(2.39) 



3. Discussion 

Figure [2] (left) shows the theoretical eigenvalues of the slip- length tensor b c g for a slip- 
ping area fraction </>2 = 0.5 as a function of the slip length b calculated with Eqs. (|2.22|) 
and (|2.35[) . In addition, we plot the data for tilted stripes (9 = 7r/4). Also included in 
Figure [5] (left) are results of a numerical solution of the Stokes' e quations performed by 
C. Co ttin-Bizonne and C. Barentin using the method developed in lCottin-Bizonne et al. 
( 2004I ). The agreement between a theory and simulation data is very good for all 4>2 
and b/L, but at b/L = 0(1) there is some small discrepancy, suggesting that our for- 
mulas slightly underestimate the effective slip, which is likely due to a simplification of 
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Eq. (|2.19p . The same trends were observed for other values of cf>2- Still, our analytical 
expressions for the effective slip, Eqs. f)2.22[) and (|2.35|) . appear to be surprisingly accu- 
rate, especially taking into account their simplicity. The same remark concerns the use 
of tensorial formula, Eq. (|1.2l) . 

Our results imply that flow past stripes is controlled by the ratio of the local slip 
length b to texture period L. At b/L ^> 1 our expressions for b e e turn to Eqs. (jl.ljl 
suggested earlier for a perfect local slip. As expected, the effective slip decreases when 
b/L = 0(1) and smaller. Interestingly, this ratio also controls the anisotropy of the flow. 
Indeed, combining Eqs. (|2 . 22[) and (|2.35[) we get 
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If b/L 3> 1, the effective slip for parallel stripes, b\ s , is twice that of perpe ndicular stripes, 
b^jf, as it was in case of a perfect s lip =00) at the liqud-gas in t erface ( Lauga fc Stond 
20031 ICottin-Bizonne et all 12004 ISbragaglia fc Prosperettl 120071: iBahga et aLll2010l) . In 



this case surface anisotropy leads to a truly tensorial effective slip. However, anisotropy of 
the flow decreases with a decrease in b/L, and at small b/L we get b^ ~ b. In other words, 
at small local slip we predict simple surface-averaged, isotropic flows (independent of 
orientation) , which means despite the fact that the local slip varies in only one direction, 
the effective slip is scalar. These unexpected results are summarized in Fig. [5] (right). 
This finding can be understood by using the following simple arguments. Following the 
advice of H.A. Stone (private communication), let us consider the average fluid velocity 
(u s ) on the SH surface. According to boundary condition ()2.3[) 
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For transverse flow this expression takes the form 
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dy 



dx 



dui 
dy 



dx. 



(3.3) 



where Cq = = const is, obviously, independent of the relative orientation of 

stripes in respect to a pressure gradient, since uq represents the solution of the problem 
for a smooth homogeneous surface. The same arguments apply in longitudinal case, where 
the only difference would be the integration over z instead of x. Therefore, when 6 is a 
small value (b/L = 0(e)), the second term in (|3.3[) may be neglected as an infinitely 
small value of higher (second) order because u\ oc e, and, thus, 



(bar ) b _>o ~ H2 + O (e 2 ) 



(3.4) 



is independent of an external force direction. The anisotropy of the effective slip is de- 
termined by the second integral term in (|3.3|) . which dominates when b/L — 0(1) and 
larger. These results suggest that both the value (upper limit) of the effective slip length 
and anisotropy of the flow are controlled by the smallest characteristic length of the 
problem (in our case, b or 8). 

Finally, we present average velocity profiles in longitudinal (9 = 0) and transverse (6 = 
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y/L y/L 

Figure 3. Average velocity profiles ((j>2 = 0.8) for a longitudinal (left) and transverse (right) 
flow. From top to bottom b/L = 1000, 1 and 0.1 

7r/2) configurations for different values of local slip length b (Fig. [3]). Mean flow remains 
two-dimensional and parabolic when the driving force is applied in main directions, yet 
both the average slip velocity at y = and the maximal velocity value at the middle of 
the channel depend on b. For arbitrary 6 the flow is essentially three-dimensional as the 
orthogonal velocity component appears due to the tensorial effective boundary condition 
Ea. (|2~36)l . 



4. Conclusion 

We have analyzed pressure-driven flow over striped SH surfaces. Unlike the previous ap- 
proach, we have obtained general analytical solutions for any value of local partial slip. We 
have confirmed that the hydrodynamic response of a striped slipping surface is generally 
anisotropic. Our main conclusion is that both effective slip and flow anisotropy are con- 
trolled by the ratio of local slip at th e gas area to texture si z e. When this ratio is la r ge, ou r 
results are closely related t o thos e of iLauea k, Stond (120031 ) : ICottin-Bizonne et al\ (|2004h : 
Sbragaglia fc Prosperettil ( 2007 ): Bahga et all < 201ph . and surface anisotropy leads to 



anisotropy of effective slip. For a small ratio we predict not only a decrease in the effec- 
tive slip, but also a different, isotropic response of the striped SH surface. 
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